^AGU.PUBLICATIONS 



Journal of Geophysical Research: Space Physics 


RESEARCH ARTICLE 

10.1 002/201 3JA01 9738 


Key Points: 

• Test particle tracing reproduces 
energetic electron events 

• Perpendicular anisotropy domi- 
nates further out, parallel anisotropy 
closer in 

• PSBL/lobe source dominates event 
via reconnection 


Correspondence to: 

J. Birn, 

jbirn@spacescience.org 


Citation: 

Birn, J., A. Runov, and M. Hesse 
(2014), Energetic electrons in dipo- 
larization events: Spatial properties 
and anisotropy, J. Geophys. Res. 
Space Physics, 119, 3604-3616, 
doi:1 0.1 002/201 3JA01 9738. 


Received 23 DEC 2013 
Accepted 8 APR 2014 
Accepted article online 3 MAY 2014 
Published online 27 MAY 2014 


Energetic electrons in dipolarization events: 

Spatial properties and anisotropy 

J. Birn 1,2 , A. Runov 3 , and M. Hesse 4 

1 Space Science Institute, Boulder, Colorado, USA, 2 Guest Scientist, Los Alamos National Laboratory, Los Alamos, New 
Mexico, USA, 3 Department of Earth, Planetary, and Space Sciences, University of California, Los Angeles, California, USA, 
4 NASA Goddard Space Flight Center, Greenbelt, Maryland, USA 


Abstract Using the electromagnetic fields of an MHD simulation of magnetotail reconnection, flow 
bursts, and dipolarization, we further investigate the acceleration of electrons to suprathermal energies. 
Particular emphasis is on spatial properties and anisotropies as functions of energy and time. The simulation 
results are compared with Time History of Events and Macroscale Interactions during Substorms 
observations. The test particle approach successfully reproduces several observed injection features 
and puts them into a context of spatial maps of the injection region(s): a dominance of perpendicular 
anisotropies farther down the tail and closer to the equatorial plane, an increasing importance of parallel 
anisotropy closer to Earth and at higher latitudes, a drop in energy fluxes at energies below ~10 keV, 
coinciding with the plasma density drop, together with increases at higher energy, a triple peak structure 
of flux increases near 0°, 90°, and 180°, and a tendency of flux increases to extend to higher energy closer 
to Earth and at lower latitudes. We identified the plasma sheet boundary layers and adjacent lobes as a 
main source region for both increased and decreased energetic electron fluxes, related to the different 
effects of adiabatic acceleration at high and low energies. The simulated anisotropies tend to exceed the 
observed ones, particularly for perpendicular fluxes at high energies. The most plausible reason is that 
the MHD simulation lacks the effects of anisotropy-driven microinstabilities and waves, which would 
reduce anisotropies. 


1. Introduction 

Energetic (i.e., suprathermal, tens to hundreds of keV) particle flux increases are a typical signature of dipo- 
larization events in the near magnetotail, associated with substorms [e.g., Parks and Winckler, 1968; Lezniak 
et al, 1 968; Arnoldy and Chan, 1 969; Baker et ai, 1 978; Belian et at., 1981] and other activity [e.g., Nakamura 
et at., 2002; Runov et at., 2009; Sergeev et at., 2009; Deng et at., 2010]. These events are characterized by a 
transient or permanent rapid change from stretched, tail like, to more dipolar magnetic fields. They may 
consist of an earthward moving brief dipolarization pulse or "dipolarization front," typically associated with 
an earthward flow burst [e.g., Nakamura et at., 2002; Runov et at., 2009; Sergeev etal., 2009] or, typically closer 
to Earth, an increase of B z toward a more persistent higher level [ Nakamura et at., 2009]. 

The probable cause of the acceleration and energetic particle flux increase is the cross-tail electric field asso- 
ciated with the dipolarization. This conclusion is supported particularly by test particle tracing in assumed 
or simulated electromagnetic fields of earthward moving electric field pulses associated with localized flow 
channels [Birn et at., 1 997, 1 998, 2004; Li et al., 1 998; Zaharia et al., 2000; Ashour-Abdalla et al., 2011; Gabrielse 
et al., 2012; Pan et al., 2014]. 

Using the three-dimensional time-dependent electric and magnetic field from MHD simulations of near-tail 
reconnection and flow bursts as basis for test particle studies, we have previously demonstrated major 
acceleration mechanisms and source regions and provided insights into spatial and temporal variations 
of electron and ion fluxes [ Birn et al., 1 997, 2004, 201 3]. In the present paper we further investigate elec- 
tron acceleration and flux properties, extending the results of Birn et al. [2013] (Paper 1) on the basis 
of the same MHD simulation [Birn et al., 2011]. The simulation results are compared to Time History of 
Events and Macroscale Interactions during Substorms (THEMIS) observations, based on a previously pub- 
lished event with a data selection discussed in section 2. In section 3 we summarize major features of 
the simulation approach and earlier findings, illustrating typical orbits and source regions of accelerated 
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electrons. Section 4 provides an extended overview of spatial characteristics of injection fronts, while 
section 5 presents pitch angle dependencies of the electron fluxes and section 6 shows energy depen- 
dencies in energy-time spectrograms and energy flux distributions, providing representations closest to 
empirical results. 

2. Observations 

In this paper, we use multipoint observations of a dipolarization front observed by THEMIS probes 
[Angelopoulos, 2008] on 27 February 2009 [Runov et al., 2009]. Spectra and energy-time diagrams were taken 
from two THEMIS probes, THC, located atx « -1 7 R E , and THD, located at x « -11 R E , between 0750 and 
0754 UT. Data from two plasma instruments, the electrostatic analyzer (ESA) [Me Fad den et al., 2008], which 
provides ion and electron distribution functions in the 5 eV to 25 keV range with A E/E < 20%, and the Solid 
State Telescope (SST) [Angelopoulos, 2008], which provides high-energy (30 keV to ~1 MeV) ion and elec- 
tron fluxes with A E/E « 30%, are combined in the spectrograms. Both instruments provide one distribution 
function per probe spin (~3 s). Magnetic field measurements are provided by the Fluxgate Magnetometer 
(FGM) [Auster et al., 2008], 

3. Simulation Approach and Earlier Findings 

In this section we summarize some findings of the earlier investigations [Birn et al., 1 997, 2004, 

2013], illustrating acceleration mechanisms and source regions. The basis of the investigation are the 
three-dimensional MHD magnetic and electric fields of a simulation of near-tail reconnection and field col- 
lapse [Birn et al., 2011], which were saved at intervals of 1 Alfven time (800 time steps, corresponding to 
~10s). The simulation included a period of external driving, which lasted until f = 61 (in dimensionless time 
units of ~10s), causing the formation of a thin embedded sheet of intensified current density in the near tail. 
After t = 61 finite resistivity was imposed, which led to the initiation of reconnection (formation of a neutral 
line) at f « 90. However, fast reconnection and the generation of fast flows and strong electric fields did not 
occur until after t « 125, when reconnection proceeded to the lobes. The speed up of the evolution appar- 
ently was caused by an interaction or reconnection with ballooning-/interchange-type modes, enabled by 
the entropy reduction in the reconnected lobe magnetic flux tubes [Birn et al., 2011]. In the following, we 
will refer to the configuration at t = 61 as the "initial state." 

A realistic conversion of the dimensionless units of the simulation can be based on a magnetic field unit of 
12.6 nT, velocity of 1000 km/s, and length unit of ~1 .5 R E , leading to a time unit of 10 s and electric field 
of 12.6 mV/m. Using these units, the inner boundary of the simulation box corresponds tox « -7.5 R E and 
the outer boundary to x « -97.5 R E . In the following, we will mostly use dimensional units but keep the 
dimensionless times, for better reference to the MHD simulation paper. 

Particle orbits were integrated numerically, using linear interpolation of the MHD fields in time. Electron 
orbits were integrated using a combination of full orbits and gyrodrift orbits, based on the conservation of 
the magnetic moment /<. The switch between the two types of orbits was based on the magnitude of an 
adiabaticity parameter, representing the square root of the ratio between field line curvature radius and 
gyoradius [Buchner and Zelenyi, 1 989; Birn et al., 2004]. Drift orbit continuity requires a cubic spline inter- 
polation in space. To avoid artificial local maxima and minima, a monotonicity-conserving algorithm was 
employed [Hyman, 1983]. 

A backward tracing technique [Curran and Goertz, 1989] was used to integrate orbits from selected "final" 
times and locations backward until they reached the initial state or a boundary of the simulation box, using 
mirroring at the inner boundary. We then applied Liouville's theorem of the conservation of phase space 
density f along a phase space trajectory to calculate f values and differential energy fluxes at the final des- 
tination from the initial and boundary values at the "source" locations. A full distribution function at the 
chosen final location and time can then be obtained by varying the final energy and pitch angle. Time pro- 
files of fluxes at selected locations, energies, and pitch angles are obtained by varying the times from which 
particles are launched backward. 

Since the MHD simulation does not provide electron information (other than density, which should equal 
ion density), we must make particular choices in defining the initial and boundary electron distributions at 
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(a) 


(b) 



t (1 Os) 


t (1 Os) 


Figure 1 . Two different characteristic orbits (heavy blue lines), indicating the two main source regions of accelerated 
electrons: (a) near-equatorial drift orbit of an electron drifting into the acceleration region from the dusk flank, overlaid 
on the cross-tail electric field (color) and contours of the magnetic field component B z (solid black lines), and (b) sections 
of a bounce orbit projected into the x,z plane, overlaid on the cross-tail electric field (color) and magnetic flux contours 
(solid black lines). Figures la and 1 b (bottom) show the electron energies as function of time. 


the source locations. For the results reported here, we imposed isotropic kappa distributions [ Vasyliunas , 
1 968; Christon et at., 1 988, 1 989] 


f(E) <x('\ + E , , ^ (1) 

V (k — 3/2)kT e ) 

with a k value of 4.5 and chose a fixed initial and boundary temperature kT e = 0.3 keV. These values were 
largely based on a comparison with the THEMIS results (section 6). 

Figures 1 a and 1 b show two typical electron orbits, illustrating the major acceleration mechanisms as well as 
the two source regions contributing to enhanced energetic electron fluxes. The colored background shows 
snapshots of the cross-tail electric field E y , indicating the earthward propagating flow burst and dipolariza- 
tion front. Black contours in Figure la represent contours of constant B z shown at intervals of ~6 nT with the 
contour nearx = -20 representing the B z = 0 line. 
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Figure 2. Spatial distribution of 5.2 keV energetic electron fluxes in the x,y and x,z planes for two different pitch angles 
(i a = 5° and a = 85°) at two different times, (top) f = 126 and (bottom) f = 130. The heavy blue contours show the 
magnitude of the cross-tail electric field at intervals of 1 .26 mV/m increasing from 1 .26 mV/m (outermost contour). The 
dashed red lines represent the B z = 0 lines in thex.y plane. 


The heavy blue line in Figure la represents the near-equatorial orbit of an electron originating from the 
dusk central plasma sheet, earthward of the reconnection site. The electron drifts into the acceleration 
region of enhanced cross-tail electric field E y (color) from the dusk, participates in the earthward collapse, 
and becomes accelerated primarily by betatron acceleration from the increase in B z along its orbit. Similar 
entry and acceleration was found also for low pitch angle orbits that consisted of multiple bounces along 
a collapsing field line. In that case the energy increase resulted from first-order Fermi acceleration of type B 
[Northrop, 1 963], providing a slingshot effect at each equatorial crossing [Birn et al., 2004, 2013], 

Figure 1 b represents an electron that originally bounces on a closed field line extending into the more 
distant tail with an apex tailward of the reconnection site. When the field line becomes reconnected, the 
electron becomes trapped in the inner portion of the field and then participates in the field collapse while 
bouncing multiple times. (Only three of the bounces are shown in Figure 1 b.) We note that during the 
crossing of the reconnection site the adiabatic assumption breaks down, and the full electron orbit was inte- 
grated, including the nonadiabatic effects of the reconnection electric field of the resistive MFID simulation, 
which is finite at the X-line. In this particular case, however, the electric field at the reconnection site (recon- 
nection rate) was still relatively small so that the direct acceleration was small. The main acceleration in that 
case then also consisted of first-order Fermi acceleration in the collapsing magnetic field. 
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Figure 3. Same as Figure 2 but for 47 keV electron energy. 


Generally, the magnetic moment and the pitch angle may change substantially during the full orbit integra- 
tion when electrons cross the equatorial plane near the reconnection site. Therefore, the final pitch angle 
is not necessarily uniquely related to the initial pitch angle, and acceleration may involve a combination 
of different mechanisms, direct acceleration as well as subsequent betatron and Fermi acceleration. Elec- 
trons captured at a later time particularly might experience more significant direct acceleration than the one 
shown in Figure 1 b. Flowever, we found that the subsequent Fermi or betatron acceleration still were the 
dominant mechanisms. 


4. Spatial Properties of the Injection Region 

An instantaneous picture of the spatial extent of injection regions is hard to obtain from observations 
without a fleet of spacecraft. The simulations, however, permit us to put local features into the context 
of large-scale maps. A first, relatively coarse equatorial map of 90°pitch angle, 180 keV proton, and elec- 
tron fluxes was given by Birn et al. [1 998], based on test particle tracing in the MFID simulation of Birn 
and Hesse [1996]. Flere we present more detailed pictures of the energetic particle flux enhancements in 
different views. 

Figures 2 and 3 show the spatial extent of the injection regions for electrons of 5.2 keV and 47 keV, respec- 
tively, for nearly parallel (5° pitch angle) and perpendicular (85° pitch angle) fluxes in thex,y and x,z planes 
at two different times, f = 1 26 (top) and t = 130 (bottom). The black contours show the magnitude of the 
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Figure 4. Simulated electron fluxes as functions of time and pitch angle for different energies as indicated at three locations along thex axis, (top) The magnetic 
field B z and (bottom) the plasma density n at the three locations. 


cross-tail electric field E y at intervals of ~ 3 mV/m. (We note an error in Paper 1: the electric field contours in 
Figure 7 corresponded to multiples of 2.5 mV/m.) Both figures show that perpendicular energetic electron 
fluxes are mainly confined to the vicinity of the equatorial plane. Only the low pitch angle fluxes extend to 
higher latitudes along the field lines, consistent with observed anisotropies [Runov et al., 2013]. At earlier 
times, farther down the tail, the near-equatorial large pitch angle fluxes clearly dominate, whereas at later 
times, closer to Earth, the field-aligned fluxes increase in relative magnitude. 

At 5.2 keV the equatorial distributions of energetic particle fluxes are quite similar for both pitch angles. 
Flowever, at 47 keV, the perpendicular flux distribution shows remarkable differences from the parallel 
fluxes, particularly in their extension toward dusk and even behind the injection front. An inspection of the 
orbits contributing to these fluxes showed that they result from electrons that drift westward at the inverse 
B z gradient associated with the dipolarization front, while losing some of their previously gained energy, 
after being accelerated in the earthward collapse. This drift may even continue on the westward edge of the 
front, leading to a deposit of accelerated particles behind the front. 

5. Pitch Angle Dependence 

The variations of electron fluxes at various energies and locations as functions of time and pitch angle are 
illustrated in Figure 4, together with the variations of B z (top) and plasma density (bottom). First, it should 
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Figure 5. Simulated parallel and perpendicular electron fluxes as functions 
of time and energy at x = -10.5, y = 0, z = 0. (top) The magnetic field B z 
and (bottom) the plasma density n at this location. 



0754 
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Figure 6. Color-coded parallel and perpendicular differential electron 
energy fluxes (in eV/s/cm 2 /sr/eV) as functions of time and energy obtained 
by THD at x « -1 1 R E . The first panel shows the electric field component 
E y , and the second panel the magnetic field component B z . 


be noted that the MHD results (top 
and bottom) show a density decrease 
coinciding with the magnetic field 
increase, as typically observed [Runov 
et al., 2011], At low, near thermal, 
energies, the energetic electron fluxes 
reflect the behavior of the MHD 
plasma density, showing a drop when 
the plasma density is reduced. 

At intermediate (5-20 keV) ener- 
gies, we see a brief enhancement of 
fluxes just prior to the dipolarization, 
followed by a drop, except for par- 
allel fluxes around 0° and 1 80°. At 
even higher energies (20-80 keV) the 
fluxes are characterized by a triple 
peak structure at 0°, 90°, and 1 80°, 
observed, for instance, by Runov et 
al. [201 2, Figure 1 0]. A similar struc- 
ture was also shown by Wu et al. 

[201 3, Figure 4], although not com- 
mented on. Perpendicular fluxes are 
more dominant farther down the 
tail, consistent with observations by 
Wu et al. [2013]. At the most distant 
location shown, the higher energy 
electron fluxes around 90° even tend 
to precede the dipolarization. 

6. Energy Dependence 

The variations of electron fluxes as 
functions of time and energy at the 
location x = -10. 5 , y = z = 0 are illus- 
trated in Figure 5, again together with 
the variations of B z (top) and plasma 
density (bottom). To improve the 
statistics, the fluxes are averaged over 
a sliding window of five data points 
between f - 0.5 and f + 0.5 (a 10 s 
window) and six pitch angles in the 
given range. For comparison, Figure 6 
shows parallel and perpendicular 
differential energy fluxes obtained 
by THC on 27 February 2009. The 
shift in the simulated fluxes closely 
resembles observed characteristics. 

It is interesting to identify the source 
regions of the electrons contribut- 
ing to the change in energy fluxes. 
Figures 7c and 7d show the source 
locations in x of the electron orbits 
at 5° and 85° final pitch angle, con- 
tributing to Figure 5, together with B z 
from the MHD simulation (Figure 7a) 
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Figure 7. Characteristics of electron fluxes at x = -10. 5, y = z = 0: 

(a) magnetic field B z obtained from the MHD simulation, (b) omnidirec- 
tional electron distributions as function of time and energy, (c) origins of 
electrons contributing to the fluxes shown in Figure 5 for 5° pitch angle, 
and (d) the same for 85° pitch angle. Black solid lines are contours of Ig J, 
indicating the region of strongest fluxes. 


and the omnidirectional distributions 
of electrons as functions of time and 
energy (Figure 7b), summed over 
1 8 pitch angles between 0 and 90° 
and averaged over three adjacent 
times, separated by 10 s. The flux 
increases at energies above ~10 keV 
reach 2-3 orders of magnitude, con- 
sistent with results of Wu et at. [2013], 
although they can vary between 
different events. 

Orange in Figures 7c and 7d indi- 
cates electrons with source locations 
in the inner tail (x > -22.5). These 
particles dominate the distributions 
before and after the dipolarization 
event. Since the electric field is very 
small during these times, they sim- 
ply follow their gradient/curvature 
drifts around Earth. Some portion 
of these electrons thus come from 
the inner boundary where their drift 
orbits intersect on the duskside. 

It is noteworthy that these origins 
also dominate during the rise of the 
fluxes just prior to f = 130 (as well 
as during the decline shortly after 
t = 143), independent of pitch 
angle and energy. The initial rise in 
energetic particle fluxes hence is car- 
ried by electrons that have drifted 
from the dusk inner tail plane sheet 
into the acceleration region, with- 
out necessarily encountering the 
reconnection site. 


Yellow and green colors indicate origins in the more distant tail (-92.5 < x < -22.5). These particles are 
originally on closed field lines which extend to the more distant tail but become trapped when their field 
lines undergo reconnection. Blue color indicates origins at the distant boundary. Since we do not consider 
particle losses or gains at the inner boundary, orbits that are integrated backward across the reconnec- 
tion site into the plasma sheet boundary or lobe region (PSBL/lobe), where field lines intersect the distant 
boundary, are reflected at the inner boundary and stopped when they reach the distant boundary. The 
actual source of these "PSBL/lobe source" electrons might be the ionosphere as well as the plasma mantle 
converging toward the plasma sheet in the more distant tail. These particles also become trapped in the 
inner tail when their field lines undergo reconnection. 

Because of the short bounce time of the electrons, both distant tail and PSBL/lobe source particles pass 
the close vicinity of the reconnection site when their field lines undergo reconnection. They hence con- 
stitute the reconnection site source considered by Pan et at. [2014] in test particle simulations based on a 
global MFID simulation. Our orbit integration switches to full orbits, conserving the pitch angle, when they 
approach the reconnection site. The particles hence experience the full nonadiabatic effects of the recon- 
nection electric field in this region. We do not have a full quantitative comparison of the direct acceleration 
and the subsequent betatron or Fermi acceleration; however, inspection of individual orbits indicates that 
the direct acceleration is less significant than the subsequent adiabatic acceleration. 
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Figure 8. Simulated energy flux distributions before and after the dipolarization at various locations indicated in the 
panels. The top right panel shows the locations as white dots in the x, z plane superposed on the energy fluxes of 5° 
pitch angle 47 keV electrons at t = 130, as presented in Figure 3. 


At energies below ~100 keV, the main contributions during the shift of the fluxes to higher energies stem 
from particles originating from the PSBL or the adjacent lobes. This is remarkable because of the low plasma 
P of 0.01 in the lobes assumed for the MHD simulation. Obviously, this result is very sensitive to the assumed 
lobe pressure and density. The fact that the shift in energy fluxes appears realistic is an indication that the 
actual lobe pressure and density are not significantly different than assumed here. These particles con- 
tribute to both the flux increases between ~10 keV and ~100 keV, and the flux decreases at energies below 
~10 keV, despite the fact that the major portion of the particles below ~10 keV have also gained substan- 
tial amounts of energy. The simple reason for this different behavior is the fact that adiabatic acceleration 
leads to significant increases in the value of the phase space distribution function in the high-energy power 
law tail, which overcomes the reduction of the source density in the lobes, whereas at lower energies the 
latter dominates. 

At energies above ~1 00 keV, Figures 7c and 7d show a combination of sources in the inner tail, primarily 
at low pitch angles, and distant tail, particularly at higher pitch angles. The full distributions thus become 
mixed by various sources. 

A more quantitative representation of the spectral changes is provided by Figure 8, showing distributions of 
parallel and perpendicular differential energy fluxes at various locations before and after the dipolarization. 
The distributions were averaged over three adjacent energy channels and finite time intervals during which 
they were found to be reasonably stable. For instance, the before and after intervals forx = -10.5,y = z = 0 
were taken between t = 120 and f = 126 before and f = 131 and t = 142 after (see Figure 5). 

For further comparison, we show TFIEMIS spectra from 27 February 2009 at two similar locations in Figure 9. 
The simulated fluxes, particularly after the dipolarization, are not as smooth as the observed ones, as 
they are evaluated for single pitch angles, whereas TFIEMIS data are integrated over a pitch angle range 
of ±15° and energy range A E/E = 0.2. Nevertheless, one can identify similarities in the typical energy 
shift and the fact that the maxima of the perpendicular flux increases tend to be larger at larger distance, 
whereas the peaks of the parallel fluxes dominate closer to Earth (and at higher latitude, consistent with 
Runovetal. [2013]). 

Differences between the simulated fluxes and the observed ones are primarily in the perpendicular fluxes 
at high energies, where the simulated ones exceed the parallel fluxes considerably. A likely reason for this is 
the absence of kinetic waves and fluctuations in the MF1D simulation, which would likely be driven by the 
simulated anisotropies and cause their reduction. 
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Figure 9. Differential energy flux distributions before and after the 
dipolarization on 27 February 2009, obtained by THC and THD, located 
at x « -17 /?£ and x ss -1 1 R f , respectively. The fluxes are averaged 
over the times indicated in the figure. 


7. Summary and Discussion 

Using a recent MHD simulation of mag- 
netotail reconnection and earthward 
flow bursts [Birn et al., 2011] as basis 
for test particle tracing, we extended 
the investigation of the acceleration of 
electrons in the electromagnetic fields 
associated with the rapid dipolarization 
of an inner magnetotail section and the 
resulting shift and increase in energy 
fluxes. Flux estimates were based on 
backward tracing from the final locations 
and Liouville mapping from the source 
regions, assuming isotropic kappa distri- 
butions as input. The MHD/test particle 
approach successfully reproduced sev- 
eral observed injection features and put 
them into a context of spatial maps of 
the injection region(s): 

(1 ) Perpendicular anisotropies dominate 
farther down the tail and closer to the 
equatorial plane. 

(2) Parallel anisotropies become increas- 
ingly more important closer to Earth 
and at higher latitudes. The vari- 
ation with distance from Earth is 
consistent with TFIEMIS observa- 
tions by Wu et al. [2013], It might 
also be consistent with observa- 
tions by Fu et al. [201 1, 2012], who 
found an association of perpendicular 
("pancake" type) anisotropy with 
growing dipolarization events (asso- 
ciated with an increase of v x ) and 
parallel ("cigar" type) anisotropy 


with decaying events (associated with a decrease in v x ), if one assumes that growing events would be 
closer to the reconnection site than decaying events. Fu et al. [2012], however, did not find such clear 
spatial association. 


Since our imposed distribution functions at the source locations of the electrons are assumed isotropic, 
anisotropies at the final location can result either from different densities at different source regions or 
injection times of particles with different pitch angles contributing to a final distribution, or from differ- 
ent histories. Although the sources are distributed in space and time, for most part we did not identify 
significant differences in the source regions or injection times, such that the history appears to be more 
important. The anisotropy effects hence indicate that betatron acceleration is more important closer to 
the reconnection site, where the increase in the magnitude of B z is more drastic, while Fermi acceleration 
becomes more important closer to Earth. 

(3) Electron fluxes drop at energies below ~10 keV, coinciding with the plasma density drop, while they 
increase at higher energy above ~1 0 keV. We identified PSBL/lobe sources, consistent with MFID, for both 
regions up to about 1 00 keV. We also found that electrons contributing to the reduction in fluxes at lower 
energies had gained energy as well as those contributing to the increases at higher energy. The explana- 
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tion for the difference in the effect is that adiabatic acceleration raises the fluxes more significantly in the 
high-energy power law tail of the distribution functions, such that it overcompensates the decrease from 
the reduction in density of the lobe source, whereas the opposite is true at lower energies. 

(4) A triple peak structure of flux increases near 0°, 90°, and 180° was found at energies above ~10 keV. 

Such structures are observed, for instance, by Runov et al. [201 2, Figure 1 0] for energies above ~1 0 keV. 

A similar structure was also shown by Wu et al. [2013], although not commented on. Right now, we don't 
have easy explanation for this result. It seems to indicate that the effects of betatron and Fermi accel- 
eration appear to vary not only as a function of distance from the reconnection site but also tend to be 
concentrated in phase space around 90° and 0°/1 80°, respectively. This might be consistent with our 
finding that the effects also vary with distance from the equatorial plane, such that perpendicular fluxes 
tend to be spatially localized near the equatorial plane while parallel fluxes extend to high latitudes, 
causing spatial variations in anisotropies observed as well [Runov et al., 201 3]. 

(5) Flux increases tend to extend to higher energy closer to Earth and at lower latitudes. Since the particle 
acceleration is a direct consequence of the properties of the electric field in the MFID simulation, these 
effects apparently are related to the increase of E y away from the reconnection site, considerably exceed- 
ing the reconnection electric field, indicating continued acceleration up to the stopping region of 

the flow bursts/dipolarization fronts, and to a decrease of the electric field strength away from the 
equatorial plane. 

Our results shed light in particular on the sources of the enhanced (and reduced) energetic particle fluxes, 
the responsible acceleration mechanisms, and the possible reasons for anisotropies. This supports the view 
that the localized v x B electric field associated with flow bursts and dipolarization fronts, as obtained 
in MFID simulations, is indeed the major mechanism affecting particle acceleration as well as causing 
anisotropies. Flowever, it is also important to note some differences, which might indicate where the 
MFID/test particle approach needs to be supplemented by additional kinetic effects. Although not obvi- 
ous from the energy-time and pitch angle-time spectrograms, the energy flux distribution functions show 
that the simulated anisotropies tend to exceed observed ones. Also, in our modeling, the perpendicular 
high-energy fluxes tend to extend to higher energies than the parallel fluxes, which is not shown in the 
observations reported here. Posssible reasons for both are differences in magnetic geometry (resulting in 
differences between curvature and gradient drifts, affecting the relative importance of Fermi versus beta- 
tron acceleration) and the fact that the MFID simulation lacks the presence of kinetic waves, expected to 
result from anisotropy-driven microinstabilities, which would reduce the anisotropies. Another feature that 
might play a role is the fact that many observed dipolarization fronts show a strong E x electric field asso- 
ciated with the kinetic structure of the current layer at the dipolarization front. This might affect energetic 
particles as well. Flowever, in the present event the E x component was significantly smaller than E y at both 
THC and THD [ Runov et al., 2009], 

Generally, we note that different observed events may show quite different signatures. For instance, the 
event presented and modeled by Pan et al. [2014] recently seems quite different from the one presented 
here and earlier by Runov et al. [2009], showing a drop in energetic particle fluxes prior to the increase, which 
is typical also for many injection events observed at geosynchronous orbit. Such differences might be due 
to differences between different events as well as to differences between different locations within the same 
event. Despite the fact that Pan et al. [2014] modeled a particular event with their global MFID simulation, 
whereas ours is more generic, several conclusions agree well between our studies. We found that the main 
source region of energetic particles during the major part of the dipolarization event is the plasma sheet 
boundary layer or the lobes. Particles from these regions enter the reconnection site and the acceleration 
region when their field lines become reconnected, as illustrated by Figure 1 b. They are hence basically con- 
sistent with the reconnection site source assumed and imposed by Pan et al. [2014]. Also, the subsequent 
adiabatic acceleration was identified as the major mechanism for enhancing the energetic particle fluxes in 
either case. 

In addition, we found that the inner plasma sheet source (not included in the Pan et al. [2014] study) is the 
main contribution during the initial rise in the energetic particle fluxes (and the decline at later times) at 
all energies and pitch angles and at low pitch angles at energies above ~100 keV. (We should note that 
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the energy thresholds separating different regimes are not to be taken as universal numbers and probably 
depend on event size and magnitude). 
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